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Abstract 


Several parameters of transonic blade- vortex interactions (BVI) are being studied and some 
ideas for noise reduction are introduced and tested using numerical simulation. The model used is 
the two-dimensional high frequency transonic small disturbance equation with regions of 
distributed vorticity (VTRAN2 code). The far-field noise signals are obtained by using the 
Kirchhoff method which extends the numerical two-dimensional near- field aerodynamic results to 
the linear acoustic three-dimensional far-field. The BVI noise mechanisms are explained and t e 
effects of vortex type and strength, and angle of attack are studied. Particularly, airfoil shape 
modifications which lead to noise reduction are investigated here. The results presented are 
expected to be helpful for better understanding of the nature of the BVI noise and better blade 

design. 
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Introduction 


In recent years, helicopters have proved to be convenient and reliable means of 
transportation, because of their ability to take off, land and maneuver in areas inaccessible to fixed 
wing aircraft. Helicopter aerodynamics are very complex, because three-dimensional, unsteady, 
transonic, viscous, aerodynamic phenomena are involved. (An excellent survey of helicopter 
aerodynamics can be found in a recent review paper by Johnson [1].) The increasing use of 
helicopters in both military and civilian applications has drawn attention to the noise that they 
generate. In civilian applications this noise can be very annoying and in military applications 
provides early warning of a helicopter's approach. Internal noise problems make helicopters 
unattractive to prospective civilian passengers and in the military often substantially exceed 
acceptable limits. Noise in the design of new or modified helicopters is today a matter of great 
concern. The helicopter noise rules (Appendix H to Federation Aviation Regulations JAR) Part 
36) contain the helicopter noise levels for which compliance must be shown for actual flight tests. 
Depending on the maximum certified weight of the helicopter, the maximum allowable noje levels 
are between 88 and 108 EPNdB (Effective Perceived Noise Lever) for flyover, 90 and 1 10 EPNdB 
for approach, and 89 and 109 for the takeoff condition. It should be reiterated that these standards 
are applicable to all type certification applications received on or after March 6, 1986 [2]. In order 
to comply with these rules, the designer must have very good estimates of the detailed noise 
characteristics of the proposed rotorcraft. 

Among the several types of helicopter noise [3], that due to Blade Vortex- Interactions (BVI) 
is one of the most important. When present, BVI leads to high levels of vibration and to intense 
noise radiation. A great deal of computational and experimental effort is being expended to 
minimize the adverse effects of BVI, so that structural fatigue and acoustic detectability can be 
reduced, and crew and passenger comfort and noise abatement can be improved. Advances in 
those areas will lead to significant increases in productivity, public acceptance, and the 
competitiveness of the U. S. rotorcraft industry. 

BVI is the aerodynamic interaction of a rotor blade with the trailing vortex system generated 
by preceding blades as shown in figure 1 (from reference [4]). In level steady-state flight the 
induced velocity of the rotor disk tends to make all the tip vortices pass under the rotor disk. 
However, during steady descending flight or maneuvers the tip vortices are forced into the rotor 
disk plane causing strong blade-vortex interactions. The occurrence of BVI is very sensitive to 
rotor design parameters and helicopter operating conditions and has a profound impact on 
helicopter performance. When it occurs BVI is loud, impulsive in character, and tends to dominate 
the other noise sources, as shown experimentally (e.g. Schmitz and Yu [5]). 

The newly developed tilt rotor aircraft also have some complicated three-dimensional BVI 
geometries that generate far-field noise [6]. The design characteristics of the tilt rotor aircraft make 
ft a suitable candidate for short-to-medium range commercial transport applications. One of the 
important issues though, is the effect of the generated noise on the community and the passengers 
Tilt rotor aircraft have additional operational options. For example, descent can be accomplished 
with high or low nacelle tilt, high or low glide slope, different airspeeds, etc. Even for a given 
nacelle tilt and flight path, lift can be distributed differently between rotors and wing. Thus, an 
optimization study of tilt rotor operations for noise reduction should be carried out. Also, tilt rotor 
aircraft have more complicated acoustic mechanisms than helicopters for various reasons. 
Previous experimental studies show that one of the dominating mechanisms during descent is 
blade-vortex interaction. 

The noise reduction and prediction capabilities are often achievable through experimental 
programs directed at understanding the noise source mechanisms and the effects of design and 
flight parameters on noise. Flight test or wind tunnel test programs can be used, but in either case 



very accurate model for BVI noise has not yet been found for practical applications. 
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some results for friture experimental BVI cases are shown. Some of those results were also 
reported in reference [32], but a more detailed discussion is given here. 


The Numerical Method (VTRAN2) 

VTRAN2 is a code [15, 16] developed for analyzing the interactions of convected regions of 
vorticity with airfoils using transonic small disturbance theory. It is based on the ADI implicit 
scheme of the LTRAN2 code [10] with the inclusion of the high frequency term as described in 
reference [11] and the addition of regions of convected vorticity using the cloud-m-cell and 
multiple branch-cut approach. The code was modified to include viscosity [12] and monotone 

switches [13]. 

The equation for the unsteady transonic small-disturbance potential <]> is: 


A<|>tt + 2B<t> x t - C(|) X x + ^y 


( 1 ) 


where 


A = 


k 2 M2 


52/3 ’ 


B = 


kM 2 

52/3’ 


C = 


1-M 2 

6 2 /3 


(y+l)M <t) x . 


<{) is the disturbance-velocity potential, M is the free-stream Mach number, 5 is the airfoil thickness 
ratio and k is the reduced frequency coc/U 0 where co is the characteristic frequency of the unsteady 
motion, c is the airfoil chordlength and U G is the free stream velocity. The quantities x, y, t, <]) in 
equation 1 have been scaled by c, c/5 1 ^, or 1 , and c5 2 ^ U Q , respectively. The exponent m in the 
nonlinear term is chosen as m = 2, which is suitable for stronger shocks as it locates the shocks 
more accurately [33]. 

For BVI the characteristic frequency of the motion depends on the lateral distance between 
the airfoil and the vortex path. For large miss distances y v the characteristic frequency w is iyy v 
and thus k=c/y v «l. Therefore, the k 2 term in A implies that the first term (i. e. high frequency 
term) in equation (1) can be dropped giving the low frequency small disturbance theory, which 
was used in the LTRAN2 code. However, for BVI cases of interest y v =c, and thus k-1 and the 
full small disturbance equation (1) (including the high frequency term) is needed. Since time us 
nondimensionalized, the choice of k for the case of BVI is irrelevant and was taken as k-1 

(co-Uq/c), meaning that one time unit corresponds to a free stream convection of one chord. 

The boundary conditions are: 

1. Upstream outer boundary condition: 

<j>x = ^)y = 0 

2. Lateral outer boundary condition: 

4>x = tj)y = 0 

3. Downstream outer boundary condition: 


as x -> - 00 


as y -> ± 00 


( 2 ) 


(3) 
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(J) x + k (j)t = 0 as x->+« 
4. Airfoil surface boundary condition (applied at y=0): 


<}> y = ^T + k ^ as 0 < x < c (5) 

where Y± defines the airfoil surface. The boundary conditions are summarized in figure 3. 

For the viscosity calculations the viscous ramp method (wedge) is used. The viscous ramp 
model simulates the shock/boundary layer interaction by placing a wedge-nosed ramp at the base ot 
the shock to obtain the reduced shock pressure rise. The surface geometry must be augmented by 
the ramp model by adding an extra viscous term in the boundary condition of equation (5). Details 
of the calculation of that viscous term can be found in reference [34]. The ramp model was derived 
for steady-state computations. However, it can be incorporated into unsteady computations in a 
quasi-steady fashion. Thus, the model is valid for low frequencies, and its use in high frequency 
problems such as BVI can only give some qualitative information about the effect of viscosity with 
almost no additional CPU time. The more complicated and CPU time-consuming lag-entrainment 
method which was also incorporated in our code was not used, because it was not superior to the 
wedge model for unsteady cases [34]. 

The classical Kutta condition is satisfied by this small disturbance formulation. We axe 
interested in cases for which the reduced frequency range is less than 4, which is the limit for the 
application of the Kutta condition [35]. 

The pressure coefficient Cp in the unsteady small-disturbance theory is. 

C p = -25 2 / 3 ( <(>x+ k <|)t) 


In addition, the wake condition 

T x + kT t = 0 

implies that a branch cut exists between the trailing edge and the downstream outer boundary, 
across which the potential jump is, A(j) = F , where F is the circulation. 

A finite vortex core is used (cloud-in-cell method) for reasons of computational stability. -The 
core has a finite square shape limited by gridlines and the vorticity is bilinearly distributed inside. 
Thus, several branch cuts (in the x-direction) are introduced. The vortex can have a free path 
(convected by the flow) or a prescribed path (miss distance yv = constant, vortex velocity-Uo)- 
Details of the theoretical formulation were given by Chang [36] and Lyrintzis [16]. 

An alternating direction implicit (ADI) method is used for the solution of the equation, where 
the high frequency term is added in the y-sweep. An approximate factorization technique with 
monotone switches [13] is used for the steady calculation, which provides a start-up solution. 
Special care is taken for the conservative differentiation along the uneven mesh. 

A (213x199) mesh is used for the calculations. The computational mesh points are clustered 
more densely near and in front of the airfoil and then are stretched exponentially from the near 
airfoil region to about 200 chords from the airfoil in x and 400 in y-direction. More mesh points 
are added in the y-direction for the more accurate evaluation of the normal derivatives on the 
Kirchhoff surface. It was shown in the past that the (213x199) mesh leads to smoother results than 
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the (213x119) mesh! The code has a high vectorization level and the CPU time for each two- 
dimensional case on a Cray-2 computer is about 4 minutes for 800 time-marching steps. 

The numerical simulation method (VTRAN2 code) introduced above is simple and efficient. 
Also, it is able to catch the main features of transonic BVI, which are the shock motions around 
airfoils. The code was shown to agree well with other, more complex approaches including Euler 
and thin-layer Navier-Stokes computations [9], Thus, we expect that our conclusions will also 
hold if more accurate Euler/Navier- Stokes predictions are used. Three-dimensionality will, of 
course, influence the results. Some of the presented results will hold for three-dimensional cases, 
but only actual three-dimensional calculations (e. g. [37]) can show that. However, at great 
distances from the airfoil the waves become very difficult to follow because of numerical diffusion 
and dispersion errors. Thus, in order to study the far-field a more accurate method is needed 
which must also include three-dimensional effects. First, we will examine the different methods 
used in acoustic studies. 


Methods of Acoustic Analysis 

Sound can be considered as the propagating part of the full solution of some fluid-flow 
problems. In a few cases the radiated sound can be found this way. However, in most cases of 
practical interest one cannot find a full numerical solution everywhere in the flow, because of 
diffusion and dispersion errors caused by increasing mesh size in the far-field. Thus, it is often 
advantageous, or even necessary, to develop ways of determining the far-field noise from near- 
field solutions. The various approaches used to find the far-field noise, as explained by George 
[38], can be classified as follows: 

1) Full Flow Field Solution : This is a calculation of the full nonlinear flow field including far-field 

waves. An example of an application of this technique for the calculation of radiating waves 
in the mid-field from transonic blade- vortex interactions is given by Baeder et al. [25, ,26]. 
However, in order to numerically resolve the details of the acoustic three-dimensional far- 
field very fine mesh should be used which makes these computations impractical for 
technically significant problems even with today's powerful supercomputers. 

2) Acoustic Analogy: This is a nonlinear near-field solution plus the Ffowcs Williams and 

Hawkings equation. This approach starts from the calculation of the nonlinear near- and 
mid-field and the far-field is found from a linear Green's function formulation evaluated in 
terms of surface and volume integrals over the retarded time transformations of the near- and 
mid-field flow and body surfaces. We would note that there are substantial difficulties in 
including the nonlinear quadrupole term (which requires second derivatives) in the volume 
integrals, especially around shock surfaces. Thus, many investigators use near-field data 
only on the blade surface, which is less accurate since the effect of shocks is not accounted 
for. For a review of the applications of acoustic analogy the reader is referred to Farassat and 
Brentner [39]. 

3) Nonlinear N ear-field Plus Kir chhoff Method: This is a calculation of the nonlinear near- and 

mid-field with the far-field solutions found from a linear Kirchhoff formulation evaluated on 
a surface surrounding the nonlinear-field. The full nonlinear equations are solved in the first 
region (near-field), usually numerically, and a surface integral of the solution over the control 
surface gives enough information for the analytical calculation in the second region (far- 
field). This method provides an adequate matching between the aerodynamic nonlinear near- 
field and the acoustic linear far-field. An advantage of the method is that the s urfac e integrals 
and the first derivatives needed can be easily evaluated from the near-field CFD data; full 
diffraction and focusing effects are included while eliminating the propagation of the reactive 
near-field. The method is simple and accurate as it accounts for the shock-related noise in the 
far-field. This technique is the basis of this research. 
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Kirchhoff's Method 


Kirchhoff s method is an innovative approach to the rotorcraft noise problem which takes 
advantage of the mathematical similarities between the aeroacoustic equations and die equations of 
classical electrodynamics. We propose to apply the considerable body of theoretical knowledge 
regarding electrodynamic field solutions to arrive at a creative and novel solution to a difficult 
problem of technological importance. 

Kirchhoff s formula was first published in 1882. Morgans [40] derived a Kirchhoff formula 
for a moving surface using a lengthy and complicated analysis. Hawkings [41] and Monno [42- 
44] rederived that formula for some special cases using the Green's function approach. Recently, 
Farassat and Myers [45] rederived the general Kirchhoff formula for a moving deformable 
piecewise smooth surface using generalized derivatives. 

Kirchhoff s formula has also been used in the theory of light diffraction and other 
electromagnetic problems (Jones [46], Stratton [47]), as well as in problems of wave propagation 
in acoustics (Pierce [48]). Hawkings [41] used this formula to predict the noise from high speed 
propellers and helicopter rotors. Korkan et al. [49, 50] calculated high speed propeller noise using 
the Kirchhoff formulation of Hawkings [51]. Morino [52, 53] used a Kirchhoff formulation for 
the development of boundary element methods. Purcel et al. [51, 52] and Baeder [53] used a 
modified Kirchhoff method including also some nonlinear effects to calculate high speed 
compressibility noise. Some examples of Kirchhoff method implementations are also given in 
reference [28]. The Kirchhoff method was also used in BVI [16, 27, 29-32] to extend the 
numerically calculated nonlinear aerodynamic BVI results to the linear acoustic far-field. 

The Green's function approach pioneered by Morino [42-44] was chosen for its simplicity 
and because the coordinate system used is suitable for BVI calculations. The method uses a 
Green's function for the linearized governing equation to derive a representation for the solution in 
terms of its values and normal and time derivatives on a closed surface S in space, which is 
assumed to include all the nonlinear flow effects and noise sources. A full three-dimensional, 
compressible formulation is used, because in this case the Green's function is simpler and because 
the method is easily extended to include spanwise variations to model three-dimensional BVI. A 
brief discussion of the Kirchhoff formulation is given in the following paragraphs; for more details 
the reader is referred to the above references. 

A Green's function for the linearized governing equation is used to derive a representation for 
the solution in terms of its values and derivatives on a closed surface S in space, which is assumed 
to include all the nonlinear flow effects and noise sources. A full three-dimensional formulation is 
used, because the Green's function is simpler in this case, and because the method can be easily 
extended to include spanwise variations to model three-dimensional BVI. 

The Green's function G for the linear unsteady potential flow is the solution of the equation: 

v 2 g - ^fu Q dt f G= s ( x_x i >y-y 1 > z z 1 ’ t_t 1 ) 

where c G is the speed of sound, U Q is the free stream velocity, subscript 1 implies the acoustic 

source, and here 8 is the well known Dirac delta function. We should note that equation (2) is the 
same as equation (1) without the nonlinear terms. The Green's function G must satisfy the 
causality condition for hyperbolic equations 
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G = = 0 for t < t ] 




The solution of equation (2) for the subsonic case is given by 

o = _5(tH+t) (9) 

47trp 

where r Q = { (x-x') 2 +P 2 [(y-y') 2 +(z-z') 2 ] } 1/2 

[r 0 -M(x-x')] 

X= c 0 P 2 
P = (l-M 2 )!/ 2 

This result is well known (see, for example, Morino [42]). If we assume that all the acoustic 
sources are enclosed by an imaginary closed rigid surface S, then we can prove (Green s theorem) 
that the pressure distribution outside this surface S is given by: 


pCxo.yo’ 2 ©) = ' 


471 


Cf 


J_ dp 
dn G 
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l_ dR ( 
Co r oP 2 * 


dr c 
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3n, 
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) + 






2 9n 


dS 1 


o 
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( 10 ) 


A 


where " 1 " denotes a point on the Kirchhoff surface, M is the free stream Mach number, c G is the 
speed of sound, subscript o denotes the transformed values using the well known Prandtl-Glauert 
transformation: 


x 0 = x, y 0 = Py> z D = Pz 

n is the outward vector normal to the surface S, and subscript x implies the evaluation at the 
retarded time tj= t-T. 

Thus, the values of the pressure and its normal derivatives on an arbitrary surface around. the 
spariwise extent of an arbitrary flow are enough to give the far-field radiation at any arbitrary 
external point. In our work we use a rectangular box coinciding with mesh points in order to 
simplify the computation. The control volume is shown in figure 4. The pressure and its 
derivatives can be numerically calculated from an aerodynamic near-field code. 

Since Kirchhoffs method assumes that linear equations hold outside this control surface S, it 
must be chosen large enough to include the region of nonlinear behavior. However, due to 
increasing mesh spacing the accuracy of the numerical solution is limited to the region immediately 
surrounding the moving blade. As a result S cannot be so large as to lose accuracy in the 
numerical solution for the mid-field. Therefore, a judicious choice of S is required for the 
effectiveness of the Kirchhoff method. A rectangular box-shaped surface (figure 4) is used for the 
calculations. The VTRAN2 code is used to calculate the solutions on the surfaces. The x-limits 
for S were varied between 0.15 and 0.50 chords and values of 0.25 chords upstream and 
downstream of the leading and trailing edges are chosen. The y limits of S for our calculations are 
varied over a range from y s = 0.25 to 4.00 chords distance from the airfoil. Higher Mach numbers 
yield higher optimum values for y s because of the stronger nonlinearities in the larger lateral extent 


8 



z'' 

s' 

/ 

/ 

X 

X 

^ of the flow region, as expected from the scaling laws of transonic flow. The similarity parameter 
K (K=(l-M 2 )/82/3 , where d is the airfoil maximum thickness) can be used for the scaling of 
transonic flow. Figure 5 shows the dependance of ys on the value of K and was derived from 
numerical experiments with our code for a (213x199) mesh. 

Strip theory approximation is used, that is, the two-dimensional VTRAN2 solution is applied 
on different segments of the blade in a stripwise manner. Blade segments ranging from two to 
sixteen in aspect ratio are used. Usually mesh limitations keep the Kirchhoff surface close enough 
to the blade where the two-dimensional strip theory solution is still valid. It can be shown that the 
wake does not need to be included in S, since the pressure is continuous across it. The distributed 
vortex, being a region where the homogeneous equation is invalid, should be included within the 
surface S. However, the effect of including or not including the vortex is comparatively small, 
particularly as it is nearly unaccelerated while it is outside S. We verified this by trying 
calculations with surfaces which included and excluded the vortex when it was near the airfoil. By 
making calculations with or without the inclusion of the tip surfaces we found [28] that they have 
only a small effect; thus they were neglected for most of the calculations. 

Finally, some effects for non-parallel blade-vortex interactions were studied in the past 
utilizing the Kirchhoff method formulation and strip theory [29]. We can simulate a vortex of 
arbitrary shape or sweep interacting with a blade using VTRAN2 results for a parallel BVI, varying 
the delay time for the different chordwise strips. This technique is limited to small intersection 

angles A,, as the three-dimensional effects become large in higher angles. However, it is expected 
to provide adequate information in small intersection angles. The same technique of delaying times 
can be used for any arbitrarily shaped vortex (i.e. circular vortex). 

Types of Unsteady Shock Motion 

Tijdeman [54] showed experimentally, using an oscillating flap, that varying airfoil surface 
boundary conditions can give three different types of unsteady shock motion: (figure 6) 

Type A shock motion, where the shock at the rear of the supersonic region merely moves back and 
forth with concurrent changes in strength. 

Type B shock motion, where the shock moves similarly to type A, but disappears temporarily 
during the unsteady motion. 

Type C shock motion, where the supersonic region disappears, but a shock wave leaves the airfoil 
and propagates forward to the far-field. 

The above three types of unsteady shock motions heavily affect the characteristics (e. g. lift, drag) 
of all unsteady transonic flows. The type of shock motion that occurs in a given situation depends 
on the flow characteristics (e. g. free stream Mach number, airfoil shape, amplitude and frequency 
of the unsteady motion). These types of shock motion can even be observed in steady airfoils with 
severe flow separation downstream of the shock waves. Their existence in BVI has been verified 
by different experiments and calculations (e. g. Tangier's experiments [55]). 

Results and Discussion 

Some mid-field calculations for BVI are performed using VTRAN2 with a refined mesh to 
follow the waves of interest. Then the Kirchhoff method is used to examine the noise at the far- 
field. The viscous shock/boundary layer interaction is used, although the effect of this term is 
weak, as shown in reference [31]. 
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We use a NACA 64A006 airfoil; the vortex strength was Q v = 0.4 (Q v is a nondimensional 

measure of the vortex strength: Civ = 2r/cUo) and the vortex miss-distance yo = -0.5 chords, for a 
fixed vortex path. The initial vortex position is x 0 = -9.51 chords and the free stream velocity is 
one (arbitrary units) so the vortex passes below the airfoil leading edge at time T = 9.51. The 
Mach numbers of 0.875, 0.854, and 0.822 correspond to shock wave motions of types A, B, and 
C, respectively, as also shown in references [9, 27]. The three different types of the unsteady 
shock motion are thus studied. For the Kirchhoff surface (figure 4) we used a span of 8 chords, 
x s = 0.25 chords and y s = 3.5, 2.5 and 1.9 chords for the three types A, B and C, respectively. y s 
is chosen from figure 5, as explained before. 

Initially, we run a couple of cases for verification purposes. Figure 7 shows the lift 
coefficient variation with instantaneous x-location of the vortex: VTRAN2 solution, NACA 

64A006 airfoil, M^O.85, a = 0°, Q v =0.4, and y o =-0.26 chords. The result of the same for 

Euler solution can be found in reference [21]. The comparison shows that the two solutions are 
very close except the same difference on the trailing edge of airfoil, which is believed to be caused 
by the finer mesh is used by VTRAN2 on the trailing edge region. We also run the case: 

NACA0012 airfoil, M^O.80, a = 0°, Q v =0.4, and y o =-0.26. The comparison with Euler and 
Navier Stokes equations and other small disturbance formulations (from reference 56) is shown in 
figure 8 and is very satisfactory. 

1. Noise Source Mechanisms of BVI 

Figure 9 shows the Cp(T) signal of the three types of unsteady shock motion in the far-field 

(point O, r=20 chords, 0=30°), as shown in figure 4, using Kirchhoff s method. The signal for 
the higher Mach number (type A) propagates upstream slower, so it takes a longer time to arrive at 
the same point. The signal consists of two disturbances (I, IT) as also shown in references [27, 
28]. The primary disturbance I is the main BVI noise and it originates at the airfoil when the 

vortex passes below the leading edge. The secondary disturbance IT corresponds to the unsteady 
shock motion and depends on the motion of the shock wave, or even the entire supersonic region 
motion induced by the vortex passage. It originates at the airfoil at a later time, (when the vortex is 
6-8 chords behind the airfoil), and depends heavily on the type of shock motion. The time delay 
of disturbance II decreases with decreasing Mach number and disappears in subcntical cases. The 
directivity of the two disturbances is very different as will be shown later. The existence of the 
second disturbance was observed computationally by George and Chang [6] and was also verified 
experimentally by Caradonna et al. [57] and Shenoy [58] and Lent et al. [59], and computationally 
by George and Lyrintzis [27, 28], Owen and Shenoy [22] and Liu et al. [26]. For example, 
Shenoy [58] performed some Schlieren experiments for a rotor and a second disturbance was seen 
to propagate in the far-field as postulated by computational results. 

From figure 9 we can see that disturbance I increases slightly with increasing Mach number. 
We should also bear in mind that the definition of C p includes division by M 2 , so the effect of the 
Mach number is stronger than it appears in the above figures. Disturbance II also exists for type A 
and B shock motions, because it is caused by the movement of the entire supersonic pocket. It 
also appears to be decreasing as we move from type C to type A. However, if we measure it from 
peak to peak (instead of just reading the max value) it still increases, but at a lower rate than 
disturbance I. 

Figure 10 shows the lift coefficient Cl (T) for a subcritical case (M=0.6) and types A, B and 
C of shock motion. We can see that the total change of Cl is higher at higher Mach numbers. 





The drag coefficient Cd will be discussed next. It is well known that in uniform subcritical 
inviscid flow Q is zero. It should be noted though, that Cd is not zero for subcritical inviscid BVI 
because of acceleration due to the vortex . When the flow becomes supercritical, the inviscid Cd 
is higher because of the formation of supersonic pockets. This was verified by running the code 
for subcritical and supercritical cases. Cd can be easily calculated in terms of the pressure 
distribution. Figure 1 1 shows Cd (T) for a subcritical case (M=0.6) and types A, B and C of 
shock motion. The change in the Cd signal is much lower than the change in the Q. However, 
there is a significant increase as we move from subcritical to supercritical cases. Also the signal is 
positive when a shock wave is always present (e. g. type A). Finally, the Cd signal seems to catch 
some part of the second disturbance whereas Cl does not. This was also detected by Liu et al. 
[26] for some BVI cases using a thin layer Navier Stokes code, but is easier to see in the case of an 
oscillating flap [17], because there the periodicity of the motion is well-defined. 

The directivity of the noise signal in the far-field is very complicated, as has been shown by 
experimental studies. Most of these experiments are also three-dimensional, which makes them 
very difficult to compare. For example, in reference 60 it was shown for a model helicopter rotor 

that the maximum signal can have a different direction if the advance ratio \l is varied. 

Directivity is studied in a vortex (or fluid) fixed coordinate system keeping the distance from 
the vortex r v constant (r v = 50 chords). The relationship between the angles 6 and 0 V in an airfoil 
fixed and a vortex fixed coordinate system is shown in figure 12. The observer is, however, fixed 
with the airfoil. Both ways of looking at the waves produced by BVI are effective in following the 
noise produced when the vortex meets the leading edge of the airfoil and, in principle, all their 
results can be transformed into each other. However, the second, fluid-fixed way is not as 
effective in following any waves which originate at the airfoil when the vortex has been convected 
a chord length or more behind the airfoil. Such waves typically are associated with unsteady shock 
motion on the airfoil, which typically originates when the vortex is 6 to 8 chord lengths behind the 
airfoil. The vortex-oriented, or fluid-fixed, viewpoint is closer to the case of a helicopter flyover. 

* However, forward-radiated disturbances in the mid-field which are emitted at times after the vortex 

passage are sometimes difficult to observe. Some further discussion on the coordinate systems can 
be found in reference [27]. 

The Cp(t) signal for different directions is plotted in figures 12 and 13 for a subcritical case 
(M=0.6) and type C shock motion, respectively (note that span = 4 chords for this case). From 
figure 12 (M=0.6) we can see that we have only one disturbance. The amplitude of that 

5 disturbance is increasing as the direction angle 0 increases from 0° to 90° (downward directivity). 

In figure 13 we can see two disturbances. Disturbance I is getting weaker as the direction angle 0 
increases from 0® to 90® (forward directivity). Disturbance II is getting stronger as the direction 
angle 0 increases from 0° to 90° (downward directivity). The two disturbances also move closer 
as the angle 0 is increased, and finally almost merge at 0 =120°. This implies a different point of 
origin. 

In figures 12 and 13 the observer is still fixed with respect to the airfoil. The observer can be 
easily made to move with the vortex, as is implied by the use of the vortex-fixed coordinate 
system. If we do that, the results for type C are shown in figure 14. The BVI signals are slightly 

modified. For example, for 0=30° (0v=5.73°) the second disturbance is larger and is closer to 
the first one, because the observer is moving closer to the airfoil after the first disturbance. (If the 
second disturbance is emitted when the vortex is 7 chordlengths downstream, then with moving 

observer 0 is increased to 76.8° and r is almost half.) Figure 15 shows frequency spectra for the 
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above case with moving observer using Fast Fourier Transform (FFT). However, these results 
do not seem to be very accurate probably due to descretization errors. Since in most experiments a 
vortex-fixed coordinate system is used, a moving observer should be used for verification of future 
experimental results. 

We propose to regard the noise due to transonic BVT as being composed of the following two 
parts: (1) linear acoustic source (dipole); (2) non-linear acoustic source (quadrupole). 

The linear acoustic source, is associated with oscillating drag and lift, Q and Cd signals as 
shown in figures 10 and 1 1, and the corresponding noise directivity presents dipole patterns in the 
same way as does linear acoustic analysis. This behavior can be observed if we look at the 
subsonic case (e. g. M=0.6) directivity in figure 12. In this case there is only disturbance I which 
has a downward directivity, since the Cl variation is much stronger than the Cd, it is expected that 
the Cl does dominate the linear signal. 

The non-linear acoustic source is related to the shock wave variations on the airfoil lower 
surface during a BVT. The noise is related to shock strength variations and the amplitude of the 
shock wave motion. Shock strength variations occur as the vortex passes below the airfoil and can 
be detected by the significant increase in the Cd signal for the supercritical cases. They generate a 
strong portion of disturbance I that dominates the signals, especially at the horizontal direction and 
has a Cd dipole directivity. Thus we can see from figure 13, which shows the directivity of the 
type C shock motion, that disturbance I now has a forward directivity. Abrupt shock wave 
motions cause variations of the entire supersonic pocket. These variations usually take place much 
later and generate disturbance II, which can also be detected from the Qj signal. Disturbance II has 
a more complicated directivity. Trying different Mach numbers and different airfoil shapes (see for 
example reference 31), we found that disturbance II becomes maximum at an an angle between 60° 
and 90°, depending on the flow parameters. 

Thus, disturbance I depends on the shock strength and the Q variations as the vortex passes 
below the airfoil. In subcritical cases only the Cl variations are relevant, and in supercritical cases 
the shock noise dominates, especially near the horizontal. Disturbance II depends on the shock 
motion that is generated when the vortex is farther downstream (6-8 chords). The two 
disturbances can have a different point of origin, which changes their relative distance, as shown in 
figure 13. Disturbance I can be reduced by weakening the shock formation on airfoil, and 
disturbance II can be reduced by stabilizing the shock wave motion on airfoil. 

3. Some Cases for Future BVI Experiments 

Some test cases for future BVT experiments are shown here. For all test cases we have a 
NACA 0012 airfoil, vortex strength Q v =0.65, a fixed vortex path, initial vortex position -9.51 
chords. The cases are: 

Case 1: M=0.70, y o =-0.26, a=0° 

Case 2: M=0.71, y o =-0.26, a=0° 

Case 3: M=0.70, y o =-0.26, a= 0.2345° (corresponding Cl without a vortex=0.04) 

Case 4: M=0.70, y o =-0.26, a= -0.2345° (corresponding Cl without a vortex=-0.04) 

For all the cases Q, Cd, and far-field directivity signals are shown in figures 16-26. From figure 
16 we can see that the slight increase for the Mach number between cases 1 and does not influence 
Cl and Cd very much. Directivity signals (figures 17-20) for both fixed and moving observer 

show that the signal is max when 0 is between 60° and 70° which corresponds to 0 v =22°-28° . 
Finally (as will also be shown in section 4c of the results) for cases 3 and 4 we can see from 


figures 21-26, that ah increase of the angle of attack (i. e. from 0 to 0.2345) decreases the resulting 
BVI noise, whereas a decrease in the angle of attack increases the resulting BVI noise. From 
figures 23-26 we can see that the directivity patterns remain almost unchanged. 


4. Noise Reduction Techniques 

Several ideas for noise reduction of transonic BVI are introduced as follows: 

a) Airfoil shape modifications: 

The transonic BVI shock motions are mainly generated on the lower airfoil surface. Thus the 
lower surface of the airfoil needs to be modified, particularly near the leading edge, because 
disturbance I is generated as the vortex passes below the airfoil leading edge and disturbance II 
depends on the shock wave motion which is influenced by changes near the leading edge. 
Initially, as a theoretical test, we used a simple cosine function (1-cosx) as the shape addedto a 
section of the airfoil (NACA64A006) that gives the continuous slope for the modified airfoil. Then 
the more practical NACA 4-digit airfoil shape is added to the airfoil. 

Figure 27 shows the modified airfoil shape in comparison with original one. The maximum 
thickness added is 0.4% and 0.8% and the total distance of the modification is 38% from leading 
edge. It should noted here that all shapes are shown three times fatter, so that the modification can 

be detected easier. The 0.8% modified airfoil was given a small angle of attack a=0.028°, so that 
the Cl stays the same, but our conclusions are also valid at a=0°. The 0.4% modified airfoil was 
run at oc=0°, since the changes in the Cj are very small in this case. Figure 27 shows the 

corresponding far-field noise signals at 0=30° and 0=90°. We can see from these figures that both 
disturbances are reduced. However, another small disturbance is generated in between, and is 
mainly due to shape irregularities near the end of the modification interval. The 0.4% modification 
produces less noise. 

Then we moved the modification interval from 38% to 26%, and used the smaller maximum 
thickness (i. e. 0.4%). Figure 28 shows the modified airfoil and the corresponding far-field noise 

signals at 0=30° and 0=90°. The 38% modification is also included in the figure for comparison 
purposes. We can see that the 26% modification interval gives somewhat better results. 

We also modified the airfoil by adding a NACA 4-digit thickness line, with the same 
maximum thickness as before (i. e. 0.4%) on 38% region of the airfoil. The disadvantage of 
adding such a shape is that the slope of the airfoil surface is no longer continuous. Figure 29 

shows the airfoil modification and the corresponding far-field noise signals at 0=30° and 0=90°, 

respectively. For this case, we used the angle of attack oc=.04® in order to get same lift with the 
unmodified airfoil. In comparison with result from cosine modification, disturbance II decreased 
substantially, which means that the shock motion was stabilized to some degree. The additional 
small disturbance that appeared in the previous cosine-shaped modification has almost disappeared. 
From figure 29 we can see that the NACA 4-digit shape produces less noise than the cosine- 
shaped modification. 

The above results showed us that a small (e. g. 0.4% maximum thickness) modification close 
to the leading edge (e. g. 25% modification section) can lower substantially the BVI noise. The 
shape of the modification is important and should be adjusted using trial and error, so that 
aerodynamic performance does not deteriorate. 


A modified airfoil SC1095rn was developed by Sikorsky Aircraft Co. to improve the 
aerodynamic characteristics of the SC 1095 airfoil. Figures 30-32 show the time history of the Q 
and Cd with the comparison between airfoils SC1095 and SC1095m for the cases of free stream 

Mach number M oo =0.6, 0.7, and 0.8 accordingly. For all cases, the vortex strength was Q v =0.4 

the vortex miss-distance yo=-0.5 chords, for a fixed vortex path, the initial vortex position is -9.51 
chords. The following things should be noted that there are two pulses found in Cd(T) and one 
pulse in Q(T), as in previous cases. The pulses for SC 1095m seemed to appear weaker than the 

ones for SC 1095, especially when the Mach number increases from M^^.6 to 0.8. As the free 
steam velocity increase, the second pulse in Cd(T) becomes stronger, which is because stronger 
shock waves are being generated. 

The corresponding far-field acoustic signal C p (T) is shown in figures 33-37 at 30° below the 
airfoil. The first pulse in C p (T), the difference between SC1095 and SC1095m airfoils remains 

about the same when the free stream Mach number increases from M^^.6 to 0.8, and the pulse 

for SC 1095 is stronger. This could be because SC 1095m airfoil enlarges the radius of leading 
edge which reduces the effects of the vortex passage. The second pulse, is stronger for SC 1095 
compared to SC1095m and the difference increases when the free stream Mach number increases 

from M^O.6 to 0.8. It is believed that this increase is related to the shock wave which is 

generated on the airfoil lower surface. The shock wave generated in SC 1095m is stronger, but it 
does not move as much, generating less noise. 

The corresponding far-field directivity of the acoustic signals showed in figures 38-43, for a 
fixed observer. For M^O.6, the first pulse reaches the maximum value at 0=90°, and there are 
no significant differences between the directivity for SC1095 and SC1095m airfoil. Starting from 
M^O.7, there are two pulses appeared in the far-field directivity acoustic signals. The first pulse 

reaches the maximum value at 0=0° and the second pulse reaches its maximum value at 0=70°, 

while the two pulses become equally significant at 0=30°, starting at M^^O.7 to 0.8. In general, 

the noise directivity patterns are the same between the two airfoil, but the pulse amplitudes for 
SC1095 are stronger than the ones for SC1095m and the differences increase when free stream 

Mach number increases from ^^,=0.6 to 0.8. As a concluding remark we can say that, although 
SC1095m airfoil was developed to improve aerodynamics, it also generated less far-field J3 VI 
noise. 

Finally, from a study of 4-digit NACA airfoil shapes we found [31] that symmetric airfoils, 
and airfoils with a point of maximum camber further upstream, seem to give less BVI noise for 
conditions producing the same lift. 

b) Splitting one vortex into two vortices: 

This idea can be realized by splitting the trailing vortex using a small airfoil slat in front of the 
main airfoil. Some exploratory experiments [61] are currently conducted with such a 
configuration. Figure 44 shows the idea behind the two-vortex model. In order to simplify the 
problem we have assumed that the distance between the vortices remains constant Figure 44 also 
shows the effect of splitting the vortex in BVI far-field noise. We varied the distance between the 
two vortices from d=.25 to .5. It is shown from figure 44 that the longer the distance between the 
splitting vortices, the more the noise can be reduced, although the noise is not reduced as much as 
with airfoil shape modification. 



c) Changing Other Parameters 

The vortex strength is a very important parameter. If we can reduce the vortex strength by one 
half by designing a better blade tip, for example, then the corresponding far-field noise will be 
reduced greatly as shown in figure 45. However, the vortex miss distance does not have a very 
strong effect (for distances 0.25 to 0.5 chords), as shown elsewhere [28]. 

When we want larger lift (thrust), the angle of attack should be increased. In such case, figure 
46 shows one example where both disturbances are reduced significantly. This explains why the 
maximum BVI noise always occurs in a certain helicopter descent region. BVI noise can be 
reduced by a proper combination of advance ratio and descent rate during approach into a noise 
sensitive area (Schmitz and Yu [5]). 


Conclusions 

An existing numerical finite difference code VTRAN2 was used to analyze noise due to 
transonic BVI. The two-dimensional unsteady transonic small disturbance equation was solved 
numerically using ADI techniques with monotone switches, including viscous effects due to 
shock-boundary interactions and the cloud-in-cell method for the simulation of the vortex. The 
Kirchhoff s method was used to extend the numerically calculated two-dimensional near-field 
aerodynamic results to the three-dimensional linear acoustic far-field. 

The effect of the three types of unsteady shock motion (A, B and C) was investigated. The 
unsteady pressure coefficients Cp(t) showed the existence of two main disturbances. The first one 
(I) is believed to be associated with the fluctuating lift coefficient (Q) (linear noise) and the shock 
strength variations generated when the vortex passes below the airfoil (nonlinear noise); 
disturbance I has a strong forward directivity. The second (IT) is believed to be associated with the 
movement of the shock wave generated at a later time (the vortex is 6-8 chordlengths behind the 
airfoil) and has a strong downward directivity (nonlinear noise). The maximum radiation occurs at 

an angle 0 between 60® and 90® below the horizontal for an airfoil-fixed coordinate system, it 
depends on both disturbances I and II and the details of the airfoil shape. 

Modified airfoil designs can be very effective in reducing the BVI noise under the same lift 
and Mach number conditions. In particular, a small (e. g. 0.4% maximum thickness) modification 
close to the leading edge (e. g. 25% modification section) can lower substantially the BVT noise. 
The overall BVI noise can be reduced by utilizing the modified airfoil SC1095m in comparison 
with airfoil SC1095, especially for the higher Mach number cases. Thus, although SC1095m 
airfoil was developed to improve aerodynamics, it also reduced the BVI noise. Experiments 
should be conducted to verify the prediction, as well as to cover evaluation of other noise source 
effects. The shape of the modification is also very important, and should be investigated further. 

It was shown that splitting the vortex into two can reduce the noise significantly. In fact, 
both disturbances are reduced. Also, the longer the distance between the two vortices the lower 
will be the noise produced. It is shown that BVI noise will be greatly reduced by reducing the 
vortex strength which can be achieved with a better blade tip design. Also, BVI noise will be 
greatly reduced by increasing the angle of the attack, which it explains why thrust can have an 
effect on BVI noise. BVI noise can be reduced by a proper combination of advance ratio and 
descent rate during approach into a noise sensitive area. 
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Figure 1. Formulation of BVI for a four blade rotor (from Schlinker and Amiet [4]). 
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Figure 2. Two-dimensional BVI. 
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Figure 5. Empirical relationship of Kirchhoff surface length y s vs transonic similarity parameter 
K. 
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Figure 6. The three types of unsteady shock motion. 
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Figure 7. Comparison of VTRAN2 with Euler solution (from [21]); 64A006, M=0.85, Q v =0.4 
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Figure 13. BVI noise directivity for type C; r v = 50 chords, span = 4 chords, fixed observer. 


































Figure 18. BVI noise directivity for NACA 0012, Clv=0.65, yo=-0.26, M=0.70, r v = 50 chords 
span = 4 chords, moving observer. 
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Figure 21. C1(T) and Cd(T) signal for NACA 0012, Clv=0.65, yo=-0.26, M=0.70, Cl (without 
vortex)= 0.04, 0.0. 
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Figure 29. Cosine and NACA 4-digit thickness curve modification of 64A006 and far-field BVI 
noise (r v =50, 0=30°, 90°); thickness: 0.4%, modification interval: 38%. 
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Figure 33. BVI noise Cp(T) signal for SC1095 and SC1095m at M^O.6 
(r=20 chords, 0=30°, span=8 chords) 
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Figure 34. BVI noise Cp(T) signal for SC1095 and SC1095m at *^=0.65 
(r=20 chords, 0=30°, span=8 chords) 
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Figure 35. BVI noise Cp(T) signal for SC 1095 and SC1095m at M^^O.7 
(r=20 chords, 0=30°, span=8 chords) 
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Figure 36. BVI noise Cp(T) signal for SC1095 and SC1095m at 14^=0.75 
(r=20 chords, 0=30°, span=8 chords) 




































Figure 42. BVI noise directivity for airfoil SC1095m at M^O.7; r v =50 chords, span=4 chords 
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Figure 45. Effect of vortex strength on far-field BVI noise (r v =50, 0=30°). 
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Figure 46. Effect of angle of attack on far-field BVI noise (r v =50, 0=30°). 



